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Abstract 

Dynamical gluon mass generation has been traditionally plagued with seagull divergences, and all 
regularization procedures proposed over the years yield finite but scheme-dependent gluon masses. 
In this work we show how such divergences can be eliminated completely by virtue of a characteristic 
identity, valid in dimensional regularization. The ability to trigger the aforementioned identity 
hinges crucially on the particular Ansatz employed for the three-gluon vertex entering into the 
Schwinger-Dyson equation governing the gluon propagator. The use of the appropriate three- 
gluon vertex brings about an additional advantage: one obtains two separate (but coupled) integral 
equations, one for the effective charge and one for the gluon mass. This system of integral equations 
has a unique solution, which unambiguously determines these two quantities. Most notably, the 
effective charge freezes in the infrared, and the gluon mass displays power-law running in the 
ultraviolet, in agreement with earlier considerations. 

PACS numbers: 12.38.Lg, 12.38.Aw 12.38.Gc 
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I. INTRODUCTION 



The dynamical generation of a non-perturbative gluon mass has been first proposed by 
Cornwall [1], and has received significant attention over the years, both from the theoretical 
and the phenomenological point of view (see, e.g., Q]). According to this picture, even though 
the gluon is massless at the level of the fundamental Lagrangian, and remains massless to 
all order in perturbation theory, the non-perturbative QCD dynamics generate an effective, 
momentum-dependent mass, without affecting the local SU(3) C invariance, which remains 
intact. The generation of such a mass has been established by studying the Schwinger- 
Dyson equations (SDEs) 3|,J4| of QCD, in a gauge-invariant framework based on the pinch 
technique (PT) l[ O, [f], Q, J and its profound correspondence with the background field 
method (BFM) <|. 

Specifically, when studying the SDE for the PT-BFM gluon propagator, A(g 2 ), one looks 
for infrared finite solutions, i.e. with A -1 (0) > (see, e.g.. (ll. TloL U, 12]). Such solutions 
may be fitted by "massive" propagators of the form A _1 (g 2 ) = q 2 + m 2 (q 2 ), where m 2 (q 2 ) 
is not "hard", but depends non-trivially on the momentum transfer q 2 . In order to ob- 
tain massive solutions gauge-invariantly, it is necessary to invoke the well-known Schwinger 
mechanism [13I . \h\ . In particular, one assumes that the strong QCD dynamics give rise to 
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20|. These 



longitudinally-coupled composite (bound-state) massless poles 
poles play a role rather like Goldstone excitations, in the sense that they preserve the form 
of the Ward identities satisfied by the Green's functions of the theory in the presence of a 
mass, but they are not associated with the breaking of any local or global symmetry. 

When the renormalizat ion-group logarithms are properly taken into account in the SDE 
analysis, one obtains, in addition, the non-perturbative generalization of a(q 2 ), the QCD 



61n( g2+4 A f (g2) 



. The 



running coupling (effective charge), of the form |l|, ]2l|, |22] a~ x {q 2 ) 
presence of m 2 (q 2 ) in the argument of the logarithm tames the Landau singularity associated 
with the perturbative (3 function, and the resulting effective charge is asymptotically free in 
the ultraviolet , "freezing" at a finite value in the infrared, namely a: _1 (0) = 61n(4m 2 (0)/A 2 ). 



As has been emphasized in the literature 23|, |24(, the generation of a gluon mass is 



intimately connected with a variety of other relat ed p henomena, and most importantly with 



the center vortex picture of confinement 



25 



271 ] . In particular, an effective low-energy 



field theory for describing the gluon mass is the gauged non-linear sigma model known as 
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"massive gauge-invariant Yang-Mills" 23]. This model admits vortex solutions, with a long- 
range pure gauge term in their potentials, which endows them with a topological quantum 
number corresponding to the center of the gauge group [Z/y for SU(N)}, and is, in turn, 
responsible for quark confinement and gluon screening 24J Specifically, center vortices of 
are assumed to form a condensate because their entropy (per unit size) 



thickness 



m 



-i 



is larger than their action. This condensation furnishes an area law to the fundamental 



28]. 



representation Wilson loop, thus confining quarks 

The general picture described above appears to be in qualitative agreement with a 
plethora of lattice simulations, where the gluon propagators (in various gauges) reach a 
finite (non-vanishing) value in the deep infrared, as would happen in the presence of a 



' mass 



This rather characteristic behavior was already observed in early studies 



30j, 



and has been fir mly established recently (in the Landau gauge) using large-volume lattices, 



for both SU(2) |3l| and SU(3) 22, U U pure Yang-Mills (no quarks included). 

It is also important to mention that a qualitatively similar situation emerges within 
the "refined" Gribov-Zwanziger formalism, presented in j^. In this latter framework the 
gluon mass is obtained through the addition of appropriate condensates to the original 



Gribov-Zwanziger action 



36. 



371 ] . Interestingly enough, one obtains a gluon mass displaying 
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381]). as well as 



power-law running, in agreement with a variety of independent studies 
the results of the present article, as explained below. 

Since the dynamical generation of an effective gluon mass is a purely non-perturbative 
effect, its technical implementation is rather intricate, and requires the harmonious synthesis 



of several ingredients [15 



id, 
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QjQ. In particular, the exact way how the Schwinger's 



mechanism will be employed is crucial for the self-consistency of the entire picture. Turns 
out that one characteristic drawback in the realization of this dynamical scenario can in fact 
be traced back to a certain subtlety in the implementation of the Schwinger's mechanism at 
the level of the relevant SDE. 

Specifically, the massless poles necessary for triggering Schwinger's mechanism and allow- 
ing the possibility of a non- vanishing A _1 (0), enter into the SDE for the gluon propagator 
through the particular Ansatz employed for the fully-dressed three-gluon vertex. Of course, 
a physically motivated Ansatz must satisfy, in addition, the correct WI, in order to preserve 



the transversality of the gluon se 
proposed over the years |l, 
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f-energy. Even though several such Ansatze have been 



121 ] . they all suffer from a typical problem: as desired, 
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A _1 (0) does not vanish; however, its value is expressed in terms of seagull integrals, i.e. 
divergent integrals of the type J k A(k) and J k k 2 A 2 (k). This fact, in turn, introduces the 
need to make sense out of these divergences, given that one is not allowed to absorb them 
into a counterterm of the type mg(A^ v )>l^, because this would compromise the gauge in- 
variance of the original Lagrangian, which at no point is to be modified. Even though a 
variety of regularizations have been proposed in the literature cited above, it is clear that the 
appearance of these divergences, the need to regularize them, and the ambiguities resulting 
in from such a regularization, are without a doubt some of the weakest theoretical points of 
this entire construction. 

In this paper we present a more refined Ansatz for the three-gluon vertex, which com- 
pletely eliminates all seagull divergences. This new Ansatz is inspired from the photon-scalar 
vertex of scalar QED, introduced by Ball and Chiu 39(. When inserted into the gluon SDE 



obtained within the PT-BFM formalism, leads to the elimination of all seagull divergences, 
by triggering a special identity, valid in dimensional regularization (DR), yielding finally a 
non- vanishing and finite value for A~ 1 (0). 

In the context of scalar QED, the identity in question, given in Eq. ( 13. 14ft . is instrumental 
in enforcing the masslessness of the photon, in the absence of any bound-state poles, i.e. 
when the Schwinger mechanism is not in operation. Specifically, the aforementioned Ansatz 
of (39!, when incorporated into the SDE for the photon, gives rise to a A _1 (0) that is 
expressed in terms of seagull contributions, which do not vanish individually, due to the 
simple fact that the charged scalars are massive already at tree level. However, the vertex 
of Ball and Chiu is such that the divergent seagull terms appear precisely in the unique 
combination that will lead to their mutual annihilation, due to the identity of Eq. (I3.14p . 

The proposed three-gluon vertex consists of two parts: (%) a part that leads to the can- 
cellation of all seagull divergences by virtue of the identity of Eq. (I3.14p . exactly as happens 
in the scalar QED case; the only difference is that now the seagull terms in question origi- 
nate from the gluonic self-interactions, i.e. they are composed by the (effectively massive) 
gluon propagator, (ii) a part that contains massless bound-state poles, thus enforcing the 
Schwinger mechanism. It is from this second part of the vertex that, after solving the 
resulting integral equation, one finally obtains a finite value for A _1 (0). 

In addition to eliminating the seagull divergences, the use of the aforementioned vertex 
brings about a further important advantage. Specifically, the SDE for the gluon propaga- 
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tor of the PT-BFM may be separated unambiguously into two distinct but coupled integral 
equations, one governing the evolution of the effective charge (running coupling), ~g 2 (q 2 ), and 
one determining the momentum-dependence of the effective gluon mass, m 2 (q 2 ). This is to 
be contrasted with the standard procedure followed in the literature, where the SDE equa- 
tion is solved for the renormalizat ion-group (RG) invariant combination d(q 2 ) = g 2 A(q 2 ), 
which is subsequently decomposed into an effective charge and a running mass according to 
d{q 2 ) = ~g 2 '{q 2 ) I \q 2 + m 2 {q 2 )), by imposing physically motivated constraints on the form of 
~g 2 {q 2 ) and m 2 (q 2 ). This procedure suffers from the obvious ambiguity of trying to extract 
two components out of a given function; instead, the new procedure, involving two individ- 
ual equations, furnishes uniquely ~g 2 (q 2 ) and m 2 (q 2 ), and it is the d(q 2 ) that is subsequently 
obtained uniquely, by combining these two quantities. 

The present article is organized as follows. In Section [Til we review the salient features of 
dynamical gauge-boson mass generation through the Schwinger mechanism, which consti- 
tutes the cornerstone of the entire approach. We explain how the aforementioned mechanism 
must be judiciously incorporated into the SDE equations of QCD, and the crucial role played 
by the three-gluon vertex. The problem of the seagull divergences, which is endemic to all 
existing approaches, is discussed, and some examples of (not fully satisfactory) attempts 
for its resolution are mentioned. In Section II I II we turn to the instructive case of scalar 
QED, and demonstrate in detail how the seagull divergences cancel out from the SDE for 
the photon propagator, by virtue of the identity of Eq. (13.141) . which is in turn triggered by 



the vertex Ansatz of 



391 ] . A counter-example of a vertex that does not trigger the identity 



is also discussed. In Section [TV] we apply the lessons of the previous section to the case 
of (quarkless) QCD. In particular, an improved Ansatz for the three-gluon vertex is con- 
structed, which incorporates the Schwinger mechanism through the appearance of massless 
poles, and, at the same time, triggers the identity of Eq. (13.141) . leading to total seagull 
annihilation. In Section [V] we obtain the system of two coupled integral equations that 
determine the momentum dependence of two RG-invariant quantities, namely ~g 2 {q 2 ) and 
m 2 (q 2 ), for the entire range of physical momenta, i.e. from the deep IR to the deep UV. 
The system is solved numerically and the obtained solutions are discussed. Most notably, 
m 2 (q 2 ) display power-law running, in agreement with various earlier considerations. Finally, 
in Section [VI] we summarize our conclusions. In addition, in three Appendices we derive in 
detail various intermediate results used throughout the article. 
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II. MASS GENERATION AND THE PROBLEM OF SEAGULL DIVERGENCES. 



The gluon propagator, A tiv (q), in covariant gauges (in particular, linear, .Retype of 
gauges, and the BFM) has the form 



A^(g) = -i 



P^)A(g 2 )+^ 



(2-1) 



Q 

where £ denotes the gauge- fixing parameter, and the transverse projector P^q) is given by 

P^(?) = ^-^r- (2-2) 

The scalar factor A(g 2 ) is given by 

A-V)=g 2 + ffl(g 2 ), (2.3) 

where U^ iu (q) = P ^{q) n(g 2 ) is the gluon self-energy. One usually defines the dimensionless 
vacuum polarization, to be denoted by II(q 2 ), as n(g 2 ) = g 2 II(g 2 ), and thus 

A-V) = g 2 [l + II(g 2 )]. (2.4) 



As Schwinger pointed out long time ago [13j, the gauge invariance of a vector field does 
not necessarily imply zero mass for the associated particle, if the current vector coupling is 
sufficiently strong. Schwinger's fundamental observation was that if (for some reason) the 
vacuum polarization of the gauge bosons acquires a pole at zero momentum transfer, then 
the vector meson becomes massive, even if the gauge symmetry forbids a mass at the level 



of the fundamental Lagrangian [14]]. Indeed, casting the self-energy in the form of (12.41) . it 
is clear that if n(g 2 ) has a pole at q 2 = with positive residue /i 2 , i.e. II(g 2 ) = f^/q 2 , then 
(in Euclidean space) 

A-V)=g 2 + /i 2 . (2.5) 

Thus, the vector meson becomes massive, A -1 (0) = /i 2 , even though it is massless in the 
absence of interactions (g — 0). 



There is no physical principle which would preclude n(g 2 ) from acquiring a pole 40]. 
Actually, the appearance of the required pole may happen for purely dynamical reasons, and, 
in particular, without the need to introduce fundamental scalar field in the Lagrangian j4l| . 
Since bound states are expected to exist in most physical systems one may suppose that, 
for sufficiently strong binding, the mass of such a bound state will be reduced to zero, thus 



generating a mass for the vector meson without interfering with gauge invariance [la . Il6 . 
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When applying the dynamical concepts described above to pure Yang-Mills theories, 
such as quarkless QCD, one assumes that, in a strongly-coupled gauge theory longitudinally 



coupled zero-mass bound-state excitations are dynamically produced [42[]. Thus, it is clear 
that a vital ingredient for this scenario is strong coupling, which can only come from the 
infrared instabilities of a non-abelian gauge theories. The aforementioned excitations are like 
dynamical Nambu-Goldstone bosons, in the sense that they are massless, composite, and 
longitudinally coupled; but, at the same time, they differ from Nambu-Goldstone bosons as 
far as their origin is concerned: they do not originate from the spontaneous breaking of any 
global symmetry. The main role of these excitations is to trigger the Schwinger mechanism, 
i.e. to provide the required pole in the gluon self-energy, and more specifically, the gauge- 
independent n(g 2 ) obtained with the PT, thus furnishing a gauge-independent dynamical 



mass for the gluons [43 ]. 



Of course, in order to obtain the full dynamics, such as, for example, the momentum- 
dependence of the dynamical mass, one must turn eventually to the SDE that governs the 
corresponding gauge-boson self-energy. The way the Schwinger mechanism is integrated 
into the SDE is through the form of the three-gluon vertex. The latter, even in the absence 
of mass generation, constitutes a central ingredient of the SDE, and plays a crucial role in 
enforcing the transversality of the gluon self-energy. Therefore, an important requirement 
for any self-consistent Ansatz used for that vertex is that it should satisfy the correct WI 
(or STI) of the PT-BFM formulation, namely 

?%^ = AT> + g )-A^(A0. (2.6) 

In addition, in order to generate a dynamical mass one must assume that the vertex contains 
dynamical poles, which will trigger the Schwinger mechanism when inserted into the SDE 
for the gluon self-energy. 

The point is that the full realization of the procedure outlined above is very subtle. 
In particular, even though the use of a three-gluon vertex containing massless poles and 
satisfying the correct WI leads indeed to a transverse and infrared finite self-energy (i.e. 
A _1 (0) 7^ 0), as expected, the actual value of A _1 (0) has always been expressed in terms of 
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divergent integrals, of the form (see, e 
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A- 1 (0) = c 1 [ A(k) + c 2 [ k 2 A 2 (k) 

Jk Jk 



(2.7) 



where (in DR) J k = fi 2e (27r)~ d f d d k, with d = 4 — e the dimension of space-time 44|. This 
is not a problem, in principle, provided that the divergent integrals appearing on the rhs 
of (12. 7p can be properly regulated and made finite, without introducing counterterms of the 
form mg(Ay V )j4^, which are forbidden by the local gauge invariance of the fundamental QCD 
Lagrangian. However, various regularization procedures introduced in the literature have 
been eventually thwarted by all sorts of additional complications of variable severity. 

The simplest regularization possibility, for example, is to employ the usual DR trick 
for eliminating quadratic divergences, namely subtract j k k~ 2 = 0. Assuming a form 
A(k) = k 2 + m 2 (k), this standard (and completely legitimate) operation, 

1 r 1 f m 2 {k) 



A{k) 



(2.8) 



<k J k k 2 + m 2 (k) Jkk 2 J k k 2 [k 2 + m 2 (k)] 

leads to a finite integral, provided m 2 (k) drops off sufficiently fast in the UV, a feature which 
is in any case expected from a dynamically generated mass. The general problem with this 
procedure, however, is the reversal of sign that it induces jl|, which eventually clashes with 
the requirement of a positive-definite A _1 (0). 

In a recent work 12| the aforementioned procedure was refined in such a way as to 
evade the sign problem. The general idea is to eliminate the perturbative tail of A(k) by 
subtracting out DR "zeros" , using the generalized formula 

\n n k 2 



k 2 



0, n 



0,1,2,... 



(2.9) 



Specifically, for large enough k 2 , A{k 2 ) goes over to its perturbative expression, to be denoted 
by A pert (/c 2 ); it has the form 

^ 2 ) = £ a -^ (2- 10 ) 



Apert ( ' 



n=0 



where the coefficient a n are known from the perturbative expansion. Then one may use ( 12. 9ft 
to regularize the rhs of (12.71) . and obtain 



167r 2 A r -4 (0) = Cl [dyy [A(y) - A pcTt (y)} + c 2 [ dy y 2 [A 2 (y) - A 
Jo Jo 



2 

pert 



(2.11) 



which is finite (and has been shown to be positive). As explained in 12J, the obvious 
ambiguity of this procedure is the choice of the point s, past which the two curves, A(y) 
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FIG. 1: Diagrams contributing to the SDE for the photon self-energy in scalar QED. 

and Ape r t(|/), are assumed to coincide (and cancel exactly against each other). Thus, the 
actual value of A~g(0) remains largely undetermined. Even though additional qualitative 
arguments may be used to restrict the allowed interval of s, thus reaching good agreement 
with recent lattice data, from the theoretical point of view it is clear that this issue is far 
from settled. 



III. SCALAR QED AND THE SEAGULL IDENTITY 

In this section we will study some of the basic issues related to the appearance and 
cancellation of seagull divergences in the context of a theory much simpler than QCD, namely 
scalar QED. Specifically, we will study the SDE governing the photon, and we will discover 
a basic identity, which, in the absence of massless poles (i.e., with the Schwinger mechanism 
"switched off") enforces the masslessness of the photon, despite the fact that individual 
seagull contributions do not vanish. In addition, we will see through an explicit detailed 
construction that the Ansatze employed for the all-order photon-scalar vertex entering into 
the SDE are crucial for the activation of this identity. 

The SDE for the photon of scalar QED is shown in Fig. HJ It is a straightforward exercise 
to demonstrate that, by virtue of the Abelian WI's satisfied by the full vertices of the 
theory, the SDE may be truncated "loop-wise" , without compromising the transversality of 
the photon, i.e., q^ )+(d2)] = ^ng* )+( * )+( * )] = . 

At the "one-loop dressed" level the SDE for the photon self-energy reads (Fig. [2]) 

n^(g) = e 2 / VfV{k)V{k + q)Y v + e 2 / T^V(k) , (3.1) 

Jk Jk 
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fc 



ft 

'vwxA/vw^ 



(di) ^ (da) 

k 

FIG. 2: The "one-loop dressed" SDE for the photon self-energy. 

where T>(k) is the fully-dressed propagator of the scalar field. Y v is the fully dressed photon- 
scalar vertex, whose tree-level expression is given by T$ = —i(2k + q) fJ- . Moreover, the bare 
quadrilinear photon-scalar vertex is given by r]S = 2ig fJ-u . The photon-scalar vertex and 



the scalar propagator T> are related by the Abelian all-order WI 

q v T u = V-\k + q)-V-\k). 



(3.2) 



It is fairly easy to demonstrate the transversality of 11^ (g), namely that q u U flu (q) = 0. To 
that end, act with q u on the two terms on the rhs of (13. ip . and notice that, by virtue of (13.21) . 
the (contracted) first term (after appropriate shifting of the integration variable, a legitimate 
operation in DR), cancels exactly against the second. Given that U^q) is transverse, it 
assumes the form U^ u (q) = Il(g 2 )P M ^(g); thus one may determine II (q 2 ) simply by taking 
the trace of both sides of (13. ip . i.e. 

-2ie 2 



n(g 2 



d-1 



(3.3) 



'k Jk 

where Eq. (13.21) was used. 

Let us compute from ( 13.3P the one-loop expression for II(g 2 ), to be denoted by H^(q 2 ). 
We have (we are using DR throughout) 



n (1 V) 



-ie 



d 



(4fc 2 - q 2 )V (k)V (k + q)-2d V (k) 



where T> (k) = (k 



m 2 ) 



' 1 . Taking the limit q —>■ 0, we find 



n (1) (o) 



-Aie 2 
d- 1 



k 2 Vl{k) -\\ V (k) 

k A Jk 



(3.4) 



(3.5) 



Of course, there is no doubt that the photon remains massless perturbatively, i.e. we must 
have that TI^ 1 ^ (0) = 0. However, the way this requirement is realized is rather subtle: the 
rhs of (13. 5p vanishes indeed, by virtue of an identity that is exact in DR, namely 

r k 2 _ d r i 

J k {k 2 -m 2 ) 2 ~ 2 J k k 2 -m 2 ' 



(3.6) 
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or, equivalent ly, 



2m 2 



(d-2) 



1 



k 2 — m 2 



(3.7) 



<k (^ 2 — tn 2 ) 2 

The relations given in ( 13.61) and ( 13. 71) can be easily verified using the standard integration 
rules of the DR [45[. Thus, the perturbative masslessness of the photon is explicitly realized 
and self-consistently enforced within the DR. Note that Eq.( l3.6p may be cast in a form that 
is particularly suggestive for the analysis that follows, namely 



dk 2 



(3.8) 



To demonstrate (13. 8p directly, i.e. without deducing it from (13. 6p . we first go to Euclidean 
space, use spherical coordinates (see IA2j) . and integrate by parts (in d dimensions). Setting 
E — y, we have (suppressing the angular integral) 



k 2 



dyy- 



dV (y) 
dy 



dyyti-VV {y). 



(3.9) 



Evidently, dropping the surface term, an operation that can be formally justified by the 
standard analytic continuation employed within the DR (see, e.g. 46J), yields immediately 
Eq. d£BD. 

We now return to the general Eq. (13.31) . In order to analyze it further we must furnish some 
information about the form of T^. Of course, any meaningful Ansatz for must satisfy 
the WI of (I3.2p . or else the transversality of U fiU (q) will be compromised from the outset. 
The form obtained by Ball and Chiu 39j, after "solving" the WI, under the additional 
requirement of not introducing kinematic singularities, is 

F " = (k + tf-k 2 + <l) ~ + A(k ' ( l)^ k + ( D- ( lK-k-q{k + ff ) M ] , (3.10) 

where A(k, q) is finite as q — > 0. Clearly the first term satisfies (I3.2p . while the part propor- 
tional to A(k, q) is identically conserved. 

It is easy to recognize that when this latter term is inserted into ( 13.31) it yields a con- 
tribution that vanishes as q — > 0, with no additional assumptions, other than the regular 
nature of A(k,q); we will therefore neglect that term in what follows. On the other hand, 
the first term of yields 



n(g 2 ) 



ie 



d 



2 2 V(k + q)-V(k) 
^ K V> ( k + q y_ k 2 



+ 2d / V(k) 



(3.11) 
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Taking the limit of Eq. fl3.lll) as q — > 0, using that 

V{k + q) - V{k) dV{k) 
(k + q) 2 - k 2 * dk 2 

we have that 



+ 0(q 2 ), (3.12) 



n(o) 



Aie 



2 



(3.13) 



d- 1 

Of course, we must have that 11(0) = 0, given that there is nothing in the dynamics that could 
possibly endow the photon with a mass; in particular, we have not employed Schwinger's 
mechanism, i.e. we have not introduced dynamical poles, and, given the form of (13.101) . 
neither kinematic ones, which might simulate the dynamical ones at the level of the SDE 
(see below). Thus, the rhs of ( 13.131) must vanish, and therefore, we must have that 



2 dV(k) d 



.Sn^—ij™- (314) 

which is the non-perturbative generalization of ( 13.81) : its demonstration proceeds exactly in 
the same way (and under the same assumptions). 

Note a crucial point: the seagull terms appearing in (13. 8p cannot be set to zero individ- 
ually, because the scalar propagator inside them is massive (already at tree- level): the only 
way to keep the photon massless, is to employ (13. 8p . which cancels them against each other. 
For example, if the term J k V(k) on the rhs were multiplied by any factor other than (d/2) 
one would be stuck with seagull divergences. 

Let us now try a different Ansatz for T M , which, due to its special form will not trigger 
Eq. (13.141) . and thus will lead to a non- vanishing (but divergent) value for 11(0). Specifically, 
consider the vertex given by 

r M = r( ) + ^[S(A; + g)-S(A;)] , (3.15) 

where £(&) is the self-energy of the scalar field, T*~ l (k) = k 2 — m 2 + E(&). Equivalently, we 
may write 

r„ = {rf - ^[(k + qf - fc 2 ]} + ^ [V~\k + q)- V-\k)] . (3.16) 

The r M in (I3.15I) - (I3.16I) satisfies again the WI of (I3.2p . and thus, as before, the transversality 
of the vacuum polarization is guaranteed. There is an important difference, however, between 
(I3.10P and (13. 16^ : the latter contains massless poles, and thus, is capable of giving rise to a 
non- vanishing A _1 (0). 
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Indeed, substituting T M of Eq. (13.161) into (13. 3p . after straightforward algebra we obtain 



I% 2 



Using that 



Aie 2 
d-1 



V(k)V(k)(k + q) 



(k-qf 



-k 2 



(t^V(k)V(k + q) 



>k Q 

we find from (13.171) (setting d = 4) 

11(0) = ie 2 



f d —^- / V(k) 



k V 2 (k) 



2 I V{k) 

k 



k 2 V 2 (k) 



(3.17) 



(3.18) 



(3.19) 



which has the general form given in (12. 7p . Evidently, the Ansatz of (13.151) does not trigger 
Eq. (13.141) . and one ends up with a non-zero A _1 (0), which, however, is expressed in terms 
of divergent seagull-type integrals. 

It is evident from the above analysis that the massless poles, indispensable as they may be 
for generating a non- vanishing A~ 1 (0), must be incorporated into the SDE with particular 
care, or else they give rise to seagull divergences. But even without this pathology, it is 
clear that the vertex of (13.151) does not constitute an optimal Ansatz. For example, if the 
(1/g 2 ) pole is considered to be of purely non-perturbative origin (as it is supposed to), it 
has vanishing perturbative expansion, and so, to all orders in perturbation theory T M = , 
which is of course not correct. 

In the next section we will see that the correct procedure is to add to the vertex of 
(13. 101) non-perturbative pole terms, in such a way as to preserve the seagull cancellation 
implemented by Eq. (13. 14ft . and, at the same time, obtain a finite A _1 (0). 



IV. FINITE GLUON MASS GENERATION 



After having fixed the ideas in the context of a simple Abelian model, we now turn to a 
pure Yang-Mills theory. In particular, we will study the SDE of the gluon propagator in the 
case of pure (quarkless) QCD, within the PT-BFM framework. As has been explained in 
detail in the recent literature jiol. I47I. 48 1. this latter formalism allows for a gauge-invariant 



truncation of the SD series, in the sense that it preserves manifestly and at every step the 
transversality of the gluon self-energy. Specifically, for the case at hand, we will consider 
only the "one-loop dressed" part of the gluon SDE that contains gluons shown in Fig. [21 
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(«l) 



FIG. 3: The "one-loop dressed" gluonic graphs contributing to the SDE for the (background) gluon 
self-energy, Tl^{q). 

leaving out (gauge- invariantly!) the "one-loop dressed" ghost contributions and all "two- 
loop dressed" diagrams. Note that the Feynman rules used to build the SD series for the 
(background) gluon self-energy, U^ u (q), are those of the BFM [9]; in particular, the external 
gluons (distinguished by the grey circles attached to them) are treated as if they were 
background gluons. The two tree-level vertices necessary for our analysis are given in Fig. HI 
as we will see in a moment, the form of these vertices is crucial for obtaining from the SDE 
precisely the right combination of terms (and with the correct relative weights) that appears 

in dm. 

In order to reduce the algebraic complexity of the problem, we drop the longitudinal 
terms from the gluon propagators inside the integrals, i.e. we set A a/3 — > —ig a pA [49J. This 
does not compromise the transversality of Ii lll ,{q) provided that we do the same on the rhs 



/ 1 . a 




gf abc [{2k + q)^g a p + 2qpg fm - 2q a g^} , 



(4.1) 




- ig 2 \J abx f xcd {g m g ufi - g^g ua + g^g a p) 

+ f adx f xbc (g^uga/3 — g^agu/3 + 9fil39ua) 

+ f acx f xbd (g^ u g a p — g^gvc 
FIG. 4: The trilinear and quadrilinear gluon vertices in the Feynman gauge of BFM. 



(4.2) 
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of the WI satisfied by T ua g, namely we have simply 



<fT vaP = [A- 1 (k + q)-A- 1 (k)]g, 



a/3 , 



(4.3) 



instead of the full WI given in (12. 6p . 

Then, it is straightforward to show that the SDE corresponding to Fig. [3] reduces to 



2 . 

q + 



2{d 



T^A(k)A(k + q)T^ + 2d 2 U{k) 

Jk 



(4.4) 



where C\ the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)]. 
r^ /3 (g, k, —k — q) is the bare three-gluon vertex in the Feynman gauge of the BFM, given 
in Fig. HI and T^p denotes its fully-dressed version. 

The function A(q) appearing on the lhs of (14.41) is the scalar part of the gluon propagator 
in the BFM, i.e. two background gluons entering; its relation to the self-energy IL^g) is 
the same as in (I2.3p . Note that A(q) is related to the standard A(q), defined in the 
gauges, by means of the powerful identity, namely A(g)[l + G(q 2 )] 2 = A(q), where G(q 2 ) is 
an auxiliary two-point function f], 5(3] whose dynamics have been studied in detail in the 



recent literature (see, e.g. [511 ]. and references therein). To further simplify the problem, 
without compromising its essential features, we will next set G{q 2 ) = 0, i.e. we effectively 
assume that, inside the quantum loops, A(q) = A(q). Thus, in what follows we will be 
dealing with a single propagator, namely A (q), and will suppress the "hats" in order to 
reduce the notation. 



A. The three-gluon vertex 

Up until this point the analysis presented in this section is completely standard within the 
PT-BFM framework. At this point enters a new ingredient, namely the judicious Ansatz for 
the three-gluon vertex which, in addition to satisfying (14. 3 j) will allow us to use the seagull 
identity (13.141) and get a non- vanishing and finite A~ 1 (0). 

To begin with, let us first write A _1 (g) in the alternative form (in Minkowski space) 

A-\q) = q 2 H'\q) - m 2 (q) , (4.5) 
The tree-level result for A _1 (g) is recovered by setting H~ 1 (q) = 1 and m 2 = 0. 
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FIG. 5: The SDE for the three-gluon vertex. All kernels are one-particle irreducible, and the 
1/q 2 pole is not kinematic but dynamical (purely non-perturbative); physically it corresponds to 
a (composite) Goldstone mode, necessary for maintaining the local gauge invariance. 



Then, an appropriate Ansatz for T va p is given by 

■f - \ (k + q) 2 H-\k + q)-eH-\k) ] ~ (0) 

lV m°/3 - (k + q) 2 - k 2 m!i mP ' [ ' 

where the term V^p contains the non-perturbative contributions due to bound-state poles 
associated with the Schwinger mechanism. Thus, V m p represents the term containing the 
1/q 2 pole on the rhs in Fig. [5j Note that we must have 

q^p = \m 2 {k) -m 2 {k + q)]g af3 , (4.7) 

in order for the Y^p of Eq. (14 .6p to satisfy (by construction) the correct WI of (14.31) . 

The Ansatz of (14.61) mimics that of Eq. (13. lUj) to the extent that the first term contains 
the right structure to produce, when inserted into the first term on the rhs of (I4.4p . the 
derivative term appearing on the lhs of (13.141) . The rhs of (13. 14j) is already there: it is the 
second term on the rhs of (14.41) . originating directly from the seagull diagram (02). 

Notice that the first term on the rhs of H4.6[) may be expanded perturbatively, whereas 
V[iap vanishes perturbatively to all orders. Qualitatively speaking, H~ 1 {q) will have the form 
(assume for simplicity a constant rh 2 ) 

H-\q) ~g 2 [ dx \n[q 2 x(x - 1) + m 2 } + 0(g 4 ) ; (4.8) 
Jo 

perturbatively (to all orders), rh 2 = 0, and one recovers the usual one-loop logarithm g 2 ln(g 2 ) 
(displaying the typical Landau pole in the IR). Thus the role of the term V^ a p is two-fold: 
(i) it tames the Landau pole inside the dimensionless perturbative logarithm, and (ii) as can 
be seen directly from Eq. (14. 5p . it can furnish an IR-finite propagator, A _1 (0) = — m 2 (0) 
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(in Minkowski space), provided, of course, that the equation governing m 2 has non-trivial 
solutions (see next section). 

An important point related to the form of is that the transverse (i.e. identically 
conserved) component associated with the first term in Eq. (14.61) has been set to zero [it 
would correspond to the term proportional to A(k, q) in Eq. (I3.10p ]. As is well known, such 
an omission, while of little importance in the IR, leads to the mishandling of overlapping di- 
vergences in the UV. This, in turn, spoils the multiplicative renormalizability of the resulting 
SDE, which must be renormalized subtractively [see also discussion following Eq. (I5.12p . 

An analogous Ansatz for the vertex V^p may be deduced following a similar philosophy. 
For example, a simple Ansatz that captures the two essential characteristics of having a 
(composite), longitudinally coupled poles, and satisfying the WI of (14.71) is 



(4.9) 



where 



V, 



r,2 



m 2 {k) — m 2 (k + q) 



9af3 



and with the transverse part V^ a/3 satisfying 



n lx V t 



0. 



(4.10) 



(4.11) 



We emphasize that, in principle, the form of V* a/3 may not be chosen at will, but must 
ideally be determined from solving the corresponding SDE for the three-gluon vertex, shown 
schematically in Figj5j Given that this task lies beyond our present powers, in what follows 
we will treat V* a/3 as being essentially undetermined [see discussion before Eq. (I4.2ip ]. 
We can write the vertex of (14.61) equivalently as 

~A- 1 (k + q)-A- 1 (k)' 



iF 



(k + q) 2 -k 2 



r 



(0) 
ua/3 



+ V 



fiafB i 



with 



where 



V 



/wt/3 



(2k + q\ 



V m p + v* a/3 , 



m 2 (k + q) — m 2 (k) 



9af3 



(4.12) 



(4.13) 



(4.14) 



(k + q) 2 - k 2 

is a residual piece, acting as an additional (non-perturbative) vertex term, 
originating from forcing the vertex to assume the form of f)4.12p . As we will see shortly, this 



The term V; 



fj.a/3 
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last way of writing Y va ^ makes the use of the basic identity of Eq. fl3.14p immediate. Thus, 
after these rearrangements, we have that the final non-perturbative effective vertex V m p 
must be transverse, 

q^V^p = 0. (4.15) 

In summary, the vertex Ansatz proposed in (14.61) [and (14.121) ] above has three important 
properties: (i) satisfies identically the Wis of ( 14.31) . a fact that ensures the transversality 
of the resulting gluon self-energy; (ii) the pole term contained in V m p makes it possible to 
have a non- vanishing A _1 (0); (iii) triggers the basic equation (I3.14p . which, in turn, allows 
one to dispose of the seagull-type terms. Thus, as we will see in the next subsection, the 
A _1 (0) obtained from the SDE is both non-vanishing and finite. 



B. The implications for the SDE 

Let us now study the effect that the three-gluon vertex of ( 14.61) has on the SDE for A(g 2 ) 
given in ( I4.4p . Substituting for the on the rhs the expression given in ( 14.121) we obtain 
after simple algebra 



with 



n( 9 ) = (7d-8)^ ' A <* + ')- A W 



(k + q) 2 -k 2 



ig 2 C K 
2{d-l) 



+ 4d 



U(q) + Uf h (q) 



A(k + q)-A(k) d 



(k + q) 2 - k 2 



A(fc) 



(4.16) 



(4.17) 



and 



n*(g) = I^A(k)A(k + q)v^ 
ik 



= f T { ^A(k)A(k + q)\V" + {V" + V*}f a P . (4.18) 
Jk 

The term in square brackets on the rhs of (I4.17P has exactly the structure needed for em- 
ploying (I3.14p . In particular, using the notation introduced in (lAip . we can write (I4.17P 



as 



n(g) = (7d - 8) q 2 R A (q) + Ad T A (q) . (4. 19) 

Note the perfect balance of relative coefficients required for the precise term T A (q) to emerge 
from the SDE. This becomes possible within the PT-BFM framework thanks to the special 
vertices shown in Fig. HI Instead, in the conventional SD formulation (e.g., in the gauges) 
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it would be very difficult to obtain the precise combination of terms needed for implementing 
(J3H. 

Perturbatively, at one loop, H(q) of (14.181) reduces to a simple and rather familiar result 
[of course, Tlfh{q) vanishes perturbatively, to all orders]. Specifically, setting A (7c) = 1/k 2 
on the rhs of (j4.17p . it is immediate to recognize that the term in square brackets vanishes, 
since it becomes proportional to the DR integral J, k~ 2 = 0, while the first term becomes 

n (9 )=-(7 rf -8) 9 ^ F ^L_. (4.20) 

This is the one-loop contribution to the gluon self-energy coming from the graphs contain- 
ing only gluons. Remember that the ghosts have been omitted without interfering with 
the transversality of the answer; their omission amounts to having in front of the leading 
logarithm the coefficient (7d — 8) instead of (7d — 6); as a result (at d = 4) the coefficient of 
the one-loop (3 function is 10Ca/48vt 2 instead of 11Ca/48vt 2 QQ. 

Let us now turn to the basic non-perturbative features of (14.161) . Since by virtue of 
(I3.14p we have that Ta(0) = 0> ft is clear that n(0) = 0. Thus, the part of the calculation 
determining 11(g) is very similar to that of scalar QED, in the sense that it keeps the gluon 
(photon) massless. On the other hand, the term IL^,(g), not present in the scalar QED study, 
makes it possible to have A _1 (0) ^ for the gluons. 

To see this explicitly, we focus on the Tlm{q) given in (14. 18ft . The integral on the rhs of 
(I4.18P receives two contributions, one from the term containing the vertex V e [given in (14.101) 
and one from the term containing the sum {V t + V T }. Let now us assume, for simplicity, 
that the (undetermined) transverse vertex V t will be such that, when added to V r [given 
in (14.141) ]. will make the contribution from {V* + V r } to become numerically subleading 
compared to that of V . For instance, V t could be such that the total contribution from 
{V* + V T } were proportional to the terms hiq 2 ) and h(q 2 ), shown to be subleading in ( 1C6I) . 
Then, keeping only V e in ( I4.18p . we obtain 

?d f 

^fh{q) = -— \ k 2 A(k)A(k + q)[m 2 (k + q)-m 2 (k)} 

q Jk 

= ~ [ m 2 (k)A(k)A(k + q) [(k + q) 2 - k 2 ] . (4.21) 

q Jk 

Then, from the Appendix [Cj Eq. ( 1C3I) . we have that Hm(0) ^ 0, which, in turn, gives rise 
to A _1 (0) 7^ 0, as announced. 
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An important consequence of this analysis is that Eq. (I4.16P can be split unambiguously 
into two parts, one that vanishes as q 2 — > and one that does not. In fact, using (14.51) on the 
lhs of (14.161) . we can assign the two types of contributions into two separate (but coupled) 
equations, i.e. 

fH~\q) = ^-J0^IL(q), (4.22) 
™ 2 (<7) = ^0[-frM ■ (4.23) 

As we will see in the next section, the first equation will determine the momentum depen- 
dence of the effective charge, and the second the running of the gluon mass. 

V. COUPLED EQUATIONS FOR EFFECTIVE CHARGE AND GLUON MASS 

In this section we will study the system of integral equations given in Eqs. (I4.22p - (I4.23I) . 
under certain simplifying assumptions. The first step in our analysis consists in rewriting 
Eqs. ( I4.22p - (l4.23p in terms of RG-invariant quantities, which will correspond to the effective 
charge and the physical gluon mass. Then, the two coupled equations will be expressed in 
terms of these two RG-invariant quantities, and will be further evaluated. We will assume 
a spectral representation for the gluon propagator [viz. Eq. (1B7|) ]. a fact that simplifies 
enormously the form of the resulting equations. Finally, we will solve the system numerically 
and study the properties of the obtained solutions. 

A. RG-invariant quantities 

It is well-known that, due to the Abelian Wis satisfied by the PT-BFM Green's functions, 
the propagator A _1 (g 2 ) absorbs all the RG logs, exactly as happens in QED with the photon 
self-energy. Specifically, let us define the renormalization constants of the gauge- coupling 
and the effective self-energy as 

A(gV) = Z A 1/2 ^ 2 )A (q 2 ), (5.1) 
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where the "0" subscript indicates bare quantities. Then, since the renormalization constants 
above satisfy the QED-like relation 



Z g = Z A 1/2 , (5.2) 

the product 

d (q 2 ) = g 2 A (q 2 ) = g 2 A{q 2 ) = d(q 2 ), (5.3) 

retains the same form before and after renormalization, i.e., it forms a RG- invariant 
(^-independent) quantity 1]. 

For asymptotically large momenta one may extract from d(q 2 ) a dimensionless quantity 
by writing, 

d(q 2 ) = (5.4) 
where ~g 2 {q 2 ) is the RG-invariant effective charge of QCD; at one-loop 

g2{q2) = = 6 In (<7 2 /A^cd)' (5 ' 5) 

where Aqcd denotes an RG-invariant mass scale of a few hundred MeV. 

The relation given in Eq. fl5.3p is true both perturbatively and non-perturbatively. In 
order to realize it non-perturbatively, let us first set 

m 2 (q 2 ) =m 2 (q 2 )H~ 1 (q 2 ), (5.6) 

where m 2 (q 2 ) is assumed to be a RG-invariant quantity, to be identified with the dynamical 
gluon mass. Then 

A(g 2 ) = 2 f iq l] 2 . , (5.7) 
q z + m z (q z ) 

and from the requirement that g 2 A(q 2 ) must be RG-invariant we have that 

9 2 H(q 2 )=g 2 (q 2 ). (5.8) 
Therefore, we finally arrive at the RG-invariant combination 

d(q 2 )=g 2 A(q 2 )=g 2 (q 2 )A(q 2 ), (5.9) 

with 

A(g 2 )= 2 , * 2V (s.io) 
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Evidently the dimensionful RG-invariant quantity d(q 2 ) is decomposed into the product of 
two individually RG-invariant quantities, the dimensionful part [g 2 +m 2 (g 2 )] corresponding 
to a massive propagator (with a running mass), and the dimensionless g 2 {q 2 ) corresponding 
to the running coupling (effective charge). 



B. The equation for the effective charge 

Even though in principle the analysis may be carried out using systematically the for- 
mulas of Appendix (TAJ) without imposing any additional constraints on A, the presence of 
the derivatives makes the numerical treatment rather cumbersome. Instead, as shown in 
Appendix ([B]), the use of the spectral representation for A results in a spectacular simplifi- 
cation. 

Specifically, assuming that A can be written as in (IB7I) . and using the expressions in 
(IBllj) . we have that 

Ad 



„, . (7d-8) » 











J ^ 





A(z)-C 



16tt 2 



ql \ ( I"- 
azz 1 — 

q z 



1/2 



(5.11) 

The equation for the effective charge, a(q 2 ) = ~g 2 (q 2 ) / Arc , will be derived from (14.22p 
after substitution of (15. lip . At this point we go to Euclidean momenta; specifically we 



set q 



-q 2 , with q\ > the positive square of a Euclidean four- vector, and define 



the Euclidean propagator as A B (g| 
follows). Then, from (14.221) we have 

H-\q 2 ) = K + bg' 



-A(—q 2 ) (we suppress the subscript "E" in what 



Az\ 1/2 



C 



(5.12) 



where h = 10Ca/487t 2 ; the discrepancy from the correct factor b = 11Ca/487t 2 , namely the 
first coefficient of the QCD one-loop /5-function, is due to the (gauge-invariant!) omission 
of the ghost loops. The (infinite) constant K is the gluon wave-function renormalization, 
introduced in order to make the equation finite, i.e. eliminate the infinite constant C. 
Note that, as is typical in this type of SDE analysis, the renormalization is carried out 
subtractively instead of multiplicatively. This is ultimately connected with the fact that 
the transverse part of the three-gluon vertex is undetermined by the gauge technique (see 
discussion in the previous section). 



22 



The constant K may be determined from (15.121) by imposing a renormalization condition 
on the function H~ 1 {q 2 ). Specifically, using the MOM-type of condition if _1 (/i 2 ) = 1, we 
have that K is given by 



K = 1 - bg 2 



" 2/4 ( 4z \ ( 4z\ 1/2 



(5.13) 



Inserting the expression for K given in (15.131) back into (15.121) . we obtain the renormalized 
equation 



H~ (q ) = 1+bg' 



5/i 2 



1-^1 A(z) 



(5.14)" 

In order to derive the equation for the effective charge a(q 2 ) = g 2 (g 2 )/47r, use the relation 
between H(q 2 ) and ~g 2 {q 2 ) given in d 5 . 8 j) . to cast (15.141) in the form 



1 



1 



9 2 (q 2 ) 3V) 



+6 



fl\ ( 4z 



4 \V2 rf/l / Az 



Az\ l/2 

.m'7J AM 

(5.15) 



Note that again, because of the mishandling of the transverse part of the three-gluon vertex, 
the rhs of ( 15.151) is not RG-invariant. The simplest way to remedy this (by hand) is to 
replace A(z) — > A(z) of Eq. (l5.10p . Thus, we arrive at 



1 



Note that 



+6 



9 2 A 



4z 



9 2 (0) sV) 



5g 2 



-b 



4z 



dz 1 + — 1 - - A(z) - 



1/2 



/ Az 



4^V /2 
1-^1 A(2) 



4^V /2 



(5.16) 
(5.17) 



C. The equation for the gluon mass 



Let us now turn to the dynamical equation governing the evolution of the mass. From 
(147231) we obtain 

m 2 (g 2 ) = 2 -bg 2 [h{q 2 ) + h(q 2 ) + h{q 2 ) + h{q 2 )] , (5.18) 

where the terms Ii(q) are given in (1C6I) . According to the discussion in (JC]), the terms h(q 2 ) 
and Ii(q 2 ) are subleading both the IR and the UV, and may be therefore safely neglected 
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to a first approximation. Then, keeping only h{q 2 ) and Is(q 2 ), we have 
2~ 



m 2 (q 2 ^ 



: bg< 



/q j roo 

dyym 2 {y)A{y) - - dyy 2 A 2 (y)[m 2 (y)]' 



(5.19) 



The next step is to rewrite this equation in terms of the RG-invariant quantities. Using 
(ElD-flEBD, we have that 



m 2 (q 2 ^ 



h 2 



M<r) F dyy&{y)[m\y)/g\y)]- 1 - / <lyy 2 ±Hy)[m 2 (y)/f(y)}' 



(5.20) 



'0 " Jq" 

Given that g 2 (y) is expected to be a much slower varying function of the momentum com- 
pared to m 2 (y), both in the UV and the IR, we will simplify the analysis by neglecting the 
derivative [g 2 {y)}' next to [m 2 (y)}'. Then, we have that 



m 2 (q 2 ^ 



t(q 2 " 



2b 
5 



rq i poo 

A(q 2 ) / dyym 2 (y){g 2 A(y)/g 2 (y)}-- / dyy 2 [m 2 (y)]'A(y){g 2 A(y)/g 2 (y)} , 

JO 1 Jq 2 

(5.21) 

which, after using (15.91) . becomes 
m 2 (q 2 ] 



t(q 2 ) 



2b 
~5 



/q j poo 

dyym 2 {y)A{y)-- dyy 2 [m 2 (y)}'A(y)A(y) 



(5.22) 



Finally, the rhs of (I5.22p is made RG-invariant by setting A(q 2 ) — > A(q 2 ) and A(y) — > d(y), 
thus obtaining 



m 2 (q 2 ^ 



2b 
~5 



A{q 2 ) I dyym 2 {y)A{y) - - / dyy 2 A 2 (y)g 2 (y)[m 2 (y))' 



(5.23) 



'0 " Jq* 

Let us now study the behavior of the solutions of (I5.23j) for asymptotically large q 2 ; in this 
limit we set A(x) — > 1/x and A(y) — > 1/y. Then, the equation reduces to 

1 



m (q ) lng^ 



iq 2 Jo 



1 f°° 

dym 2 {y) - - dyg 2 (y)[m 2 (y)}' 

1 Jq 2 



(5.24) 



It is relatively straightforward to establish that the asymptotic solutions of (I5.24p display 
power-law running. Indeed, substituting on both sides of (I5.24p a m 2 (q 2 ) of the form 



m 2 {q 2 ) = ^{W) 



(5.25) 



it is easy to recognize that the second term on the rhs of f!5.25j) is subleading. Indeed, in 
the absence of ~g 2 {y) = (Mny) -1 the integrand is a total derivative, which yields to the 
rhs simply a term \m 2 (q 2 ); this is suppressed, because it is not multiplied by a lng 2 . The 
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presence of g 2 (y) suppresses this integral even further. Specifically, integration by parts and 
use of the second equation in (IC7I) yields 



dyg 2 (y)[m 2 (y)]' = -m 2 (x)g 2 (x) + O (1/lnx) 

which is indeed further suppressed by an extra logarithm. 
Thus, using the elementary integral 

dy 1 



/ 



y(\ny) 1+a a(\ny) a 
[first equation in (1C7j) ]. we have that (I5.25j) is a solution of (I5.24p provided that 



7= 5' 

and so, the asymptotic solution has power-law running, given by 



Ar 



™V) = -i(lnq 



2\-3/5 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



Finally, if we were to assume the approximate validity of (15.291) for the entire range of 
momenta, we can set [m 2 (y)}' ~ —m 2 (y)/y; that way, we convert (15.241) from an integro- 
differential equation to the simpler integral equation 

m 2 (q 2 ) 



g 2 (q 2 ) 



2b 
~5 



/q poo 
dyym 2 {y)A{y) + - J dyyK' 



(y)g (y)m (y) 



(5.30) 



D. Solving the system numerically 

We will next discuss the numerical solutions for the system of integral equations, namely 
( I5.16P and (I5.30p coupled together. 

We solve numerically the two coupled integral equations, renormalizing them at three 
different points, namely /i = {4, 10, 91} GeV, with a{^ 2 ) = g 2 {fi)/47T = {0.341, 0.229, 0.127}, 
respectively. In Fig. [6l we show the results for a(q 2 ); there we see clearly that the three 
curves merge practically into a single one, thus confirming numerically the /i-independence 
of <y{q 2 ), expected on formal grounds. These three curves may be accurately fitted by the 
physically motivated functional form [l|, namely 

a(q 2 ) = , (5.31) 

47rMn[(g 2 + tm2)/A2] 
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FIG. 6: Numerical solutions for the effective charge obtained from Eq. (|5.16p . renormalized at 
three different points: \i = 4GeV and a(p 2 ) = 0.341 (black curve), \x = lOGeV and a(/u 2 ) = 0.229 
(red curve), p = 91 GeV and a(p 2 ) = 0.127 (blue curve). The three curves practically coincide, 
showing that indeed a(q 2 ) is independent of the renormalization point chosen. The dashed curve 
(magenta) is the perturbative one-loop behavior, and the brown line with circles depicts the fit of 
Eq. ||EH]> , for t = 3.7. 

with t = 3.7 and A = 645 MeV [see caption of Fig. [BJ. 

In Fig. we show the dynamical gluon mass, m 2 (q 2 ), obtained as solution of Eq. (I5.30p 
at the same renormalization points of Fig. [61 Once again, this figure shows us that m 2 (q 2 ) is 
also a RG-invariant quantity, since the three curves, obtained using the three different (and 
quite disparate) renormalization points, are practically on top of each other. The behavior 
of m 2 (q 2 ) in the entire range of momenta can be accurately described by the following 
parametrization 



2 / , ,2 \ _ ln / q 2 + f(q 2 ,m 2 Q ) \ / ln ff(o,m 2 ) 



2, 

m (q 



q 2 + ml 



A 2 J / V A 2 



-3/5 



(5.32) 



where the function 



m 



f(g » m o) = Pl m Q + P2^~ 2 ' ( 5 - 33 ) 

q z + rriQ 

with pi = —1/2, p2 = 5/2, and tuq = 612 MeV. Notice that in the UV asymptotic limit the 
above expression goes over to that of Eq. (15.291) . as it should. 

Finally, we turn to the RG-invariant quantity d(q 2 ) } which appears in a natural way in 
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FIG. 7: The effective gluon mass, m 2 (q 2 ), for the same renormalization points used in Fig. ([6]). 
Evidently, the three curves merge into a single one, showing that the numerical solutions are 
independent of the renormalization point. The continuous line in magenta is the fit of Eq. (|5.32p 
with p\ = —1/2 and pi = 5/2. 



all physical process involving gluon exchange. With the help of Eq. (15. 9p we can construct 
d(q 2 ) out of the numerical solutions for a(q 2 ) and m 2 (q 2 ); the result is shown in Fig. [HJ 
Obviously, since d(q 2 ) is built out of two quantities that are individually independent of [/,, 
it too turns out to be //-independent; this property is clearly observed in Fig. [HJ 
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H = 91 GeV a(n 2 ) = 0.127 

H= 10 GeVa(n 2 ) = 0.229 

n = 4 GeV a(n 2 ) = 0.341 




1000 10000 
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FIG. 8: The renormalization-group invariant product d(q 2 ) obtained by combining the results for 
a(q 2 ) and m 2 (q 2 ) according to Eq. (|5.9p . 
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VI. DISCUSSION AND CONCLUSIONS 



In this article we have demonstrated how to obtain a finite gluon mass from the SDE of 
QCD, formulated in the PT-BFM framework. Obtaining a finite mass without the appear- 
ance of seagull divergences has been a long-standing problem, that has afflicted all related 
studies for a number of years. The key observation that leads to the solution of this problem 
is that a judicious Ansatz for the three- gluon vertex eliminates all seagull divergences by 
means of a basic identity, valid in dimensional regularization. In retrospect one realizes 
that the problem of seagull divergences is not intrinsic to this approach, but has rather 
been caused by the inadvertent mismatch of two field theoretic mechanisms, induced by an 
imperfect Ansatz for the vertex. Specifically, the Schwinger mechanism, which requires the 
appearance of massless poles in the three-gluon vertex, distorts the mechanism responsible 
for the cancellation of the seagull divergences, unless the poles enter into the gluon vertex 
in a very particular way. 

The procedure described in the present work furnishes eventually two separate but cou- 
pled equations for the QCD effective charge and the gluon mass, which, when solved simul- 

i quantities. This is a considerable improvement 



11]) where one had only one dynamical equation, 



taneously, yield a unique answer for bot 
over the existing approaches (e.g. 
determining d(q), which was subsequently decomposed according to Eq. (I5.9p . in order to 
obtain (not without a certain ambiguity) the effective charge and gluon mass. It should be 
emphasized, however, that even though the values obtained for a(0) and mo are very rea- 
sonable, they are not directly comparable with the values obtained from phenomenological 
studies, due to the fact that the gauge used in this work ("stagnant" or "generalized Feyn- 
man" ) is not the canonical Feynman gauge of the BFM, which, as is well-known, furnishes 
the PT effective charge and gluon mass. This may account, in retrospect, for the slightly 
elevated value of a(0) ~ 1 obtained here, compared to a value of 0.5 — 0.7, found in recent 
theoretical analysis [52], and various phenomenological studies 53l. 1541] 

The fully dressed three-gluon vertex T m p used in Eq. (14. 6 j) satisfies (by construction) the 
WI given in Eq. (14.31) . or, after restoring the longitudinal terms, the full WI of Eq. 



which is crucial for ensuring the transversality of the gluon self-energy. However, this Ansatz 
used for V m p is still incomplete from the point of view of full Bose symmetry, in the sense 
that it does not satisfy the correct STI when contracted with respect to the two other legs; 
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we remind the reader that the other two legs (internal lines, irrigated by virtual momenta) 
correspond to "quantum" gluons, as opposed to the external (background) gluon, where the 
physical momentum q enters. The corresponding STI satisfied when contracting with 
respect to the quantum legs is a variant of the well-known STI satisfied by the conventional 
vertex 55], and has been derived in j^]. The construction of a three-gluon vertex satisfying 



the correct WI and STIs is currently under investigation, and we hope to report the results 
in the near future. 

An additional technical issue related to the form of r m p is the omission of the identically 
conserved part, as already mentioned in the cor resp onding sections. It would be most 



interesting to extend the QED construction of 56|, |57| to the case of the three-gluon vertex 
T^p. Such a task, however, appears to be of formidable logistic complexity, given that 
there are thirteen linearly independent tensorial structures, with their corresponding form- 
factors. The techniques and the special tensorial bases introduced in 58j] may prove useful 
for simplifying such a task. 

In a similar spirit, the part V^" 1 of the vertex containing the massless poles, thus trigger- 
ing the Schwinger mechanism, should also be appropriately extended, to satisfy the correct 
WI with respect to all three legs. A prime candidate for this role is the vertex proposed 
in 5^1, given by 

V^(k h k 2 ,k 3 ) = k ^^^ P;{hW{h)- k l [m 2 (k 2 ) - m^h)} P-MF^kHcp, 

(6.1) 

where "cp" denotes "cyclic permutations". Note that a vertex such as (16. ip will furnish a 
concrete (but still not unique) expression for V t , which, in turn, will allow one to scrutinize 
some of the assumptions made in Section IIVI Of course, for self-consistency, one should 
perform the analysis in, e.g., the Landau or conventional Feynman gauges, rather then the 
"generalized" Feynman gauge employed here; calculations in this direction are already in 
progress. 
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APPENDIX A: SOME USEFUL RELATIONS 



Let us define the following quantities, 

f (k + q) - f(k) 
(k + q) 2 - k 2 



R f (q) = [ 
Jk 

T f (q) = [ k 

Jk 



2 /(* + <?)-/(*) , d 



+ o / /(*) > ( Al ) 



(k + q) 2 -k 2 2 Jk 
for an arbitrary function f(x) that is finite at the origin. 

Let us define q 2 = x, k 2 = y, (k + q) 2 = z, and let us write the (Euclidean) integration 
measure [d d k = id d k E ] in spherical coordinates 

d d k E = 2n d6(sm9) d - 2 / dy y. (A2) 
Jo Jo 



We then have that z = y + x + 2^/xycosO, and we define w = z — y = x + 2^/xy cos 9. 
Finally, recall the elementary integral 

{W m + 1 n+l \ 
\ 2 111 2 > n — Ob 
r (^) ' (A3) 
, n = 2k + 1 

Rf{q) and Tf(q) may be expanded systematically as a power series in q 2 . To that end we 
consider the Taylor expansion of f(z) around w = 0, which gives (we are assuming finite 
derivatives at the origin), 

= m + + f r\v) + ... <A4) 

w 2! 3! 

where the primes denote differentiations with respect to y. Then, one must collect pieces of 
a given order in q 2 from the various powers of w, using flA3j) . 

It is clear, for example, that when the term f'(y) on the rhs of (IA40 is inserted into Tf(q) 
generates the seagull identity (13.141) . while all remaining terms proportional to positive 
powers of w; so, 

I>(0) = 0. (A5) 

As a second example, we determine the term of Rf(q) linear in q 2 , to be denoted by Rj(q); 
to accomplish this one must collect the appropriate contributions coming from both the 
second and the third term on the rhs of ( IA4I) . Using ( 13. 14ft one then obtains 



Rf\q)=f I 
Jk 



\f"(k 2 ) + lk 2 f"(k 2 ) 
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(A6) 



or after partial integration, assuming that [y/']o° = and [2/ 2 /"]o° = (valid when /(0) is 
finite, and f(y) ~ y (or faster) at infinity), we have 

Rf{q) = C -q 2 f{U), (A7) 

where c = z/167r 2 . 



APPENDIX B: THE SPECTRAL REPRESENTATION 



Let us consider a simple massive tree-level propagator, 

1 



d m (q) 



q 2 — m 2 



(Bl) 



and set / = d m directly into ( 1A1I) . Turns out that both Rd m (q) and Td m (q) can be calculated 
exactly; specifically, using that 

d m {k + q) - d m (k) 



(k + q) 2 - k 2 



-d m (k)d m (k + q) 



it is elementary to show that 

Rd m {q) = c 



dx In ( 1 + 



q 2 x(x — 1) 



Cn 



where the cutoff-dependent constant C m is given in dimensional regularization by 



C,, 



— 7 — In 



m 



47T/i S 



(B2) 



(B3) 



(B4) 



Evidently, Rd m (0) = —cC m . As an additional check, note that the term of Rd m (q) linear in 
q 2 , obtained by Taylor-expanding flB3p . is given by 

,2 



c q 



6 m 2 



(B5) 



which coincides with the result obtained when substituting /(0) = d m (0) = in the 

general formula of (jA~l) . 

For Td m (q) we have, using the identities fl3.6p - fl3.7l) . 

q 2 + 2q ■ k 



T dm (q) = m 



cm 



k (k 2 — m?) 2 [{k + q) 2 — m 2 ] 
1 x{2x-l) 
q 2 x(l — x) + m 2 
q 2 x(x — 1 



—cm q I dx 
■I 



o 



dx In 1 + 



o 



(B6) 
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Evidently, 7d m (0) = 0, in agreement with (1A5|) . 

The results of (IB3l) - flB6l) may be used in a more gen eral way. Specifically, if we assume a 
spectral representation for the gluon propagator jlilol]]. namely 

P(A 2 ) . 



then from (1B3[) (with m — > A) we have 



dX 2 



q 2 — A 2 + ie 



(B7) 



-Ra(<?) = c 



dX 2 p(X 2 



dx In ( 1 + 



T A (g) = c dX p(X )X / da; In 1 + 



q 2 x(x — 1) 



where 



C 



dX 2 p(X 2 )C x . 



(B9) 



Note that after the renormalization of the corresponding SDE (see Section |V]) C will even- 



tually drop out. Then, the use of the following identities 

q 2 x(x — 1 



jdX 2 p{X 2 ) jf dx\xi(l + 
jdX 2 p(A 2 )A 2 jf dxlnfl + 



fill 

9 2 /4 



A 2 

g 2 x(x — r 

A 2 



\ /-« 2 / 4 / 4z\ 1/2 



</ 2 /4 / 4 ^\l/2 







9' 



1 



(BIO) 



allows one to cast Ra(q) and T&(q) again as an integral containing the gluon propagator A, 
namely 





rq 2 /* / 


#a(<?) = c 












Ta(?) = c / 


dzzll 


J 





Az\ 1/2 



(Bll) 



Note that the simple change of variables t = 4z/g 2 allows one to cast Ra(q) and T^{q) in 
the alternative form 

•i 



-Ra(<?) = c 



(g 2 /4) / dt(l-t) 1/2 A(tq 2 /4) - C 
Jo 

T A (q) = c{q 2 /Af f dtt(l-ty/ 2 A(tq 2 /4), 
Jo 



(B12) 



which makes the identification of the IR behavior of these quantities immediate, and is 
particularly useful for their numerical treatment. 
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APPENDIX C: THE MASS EQUATION 



Let us consider (in Euclidean space) the integral appearing on the rhs of (14.211) . to be 
denoted by I(q)- We have 

I(q) = -J m 2 (k)A(k)A(k + q)[(k + q) 2 -k 2 ], (CI) 

E 



1 Jk 



(C2) 



which, with the notation introduced in Appendix (jA]), reads 

I(x) = [ m 2 (y)A(y)A(z) l + ^2cosfl 

Jk E L V X 

Then, expand A(z) = A(y) + wA'(y) + and collect the terms that survive the angular 
integration, to obtain 

1(0) = j m 2 (k 2 )A(k 2 )[A(k 2 )+Ak 2 A'(k 2 )cos 2 e] 

J k E 

= j m 2 (k 2 )A(k 2 )[A(k 2 ) + k 2 A'(k 2 )] 

J k E 

= -U k 2 A 2 (k 2 )[m 2 (k 2 )}'. (C3) 

Note that a monotonically decreasing mass, [m 2 (k 2 )]' < 0, guarantees that 1(0) > 0, or, 
equivalently, the positivity of m 2 (0) in Euclidean space. 

To write I(x) in a form suitable for solving the corresponding dynamical equation, first 
split the radial integration into two intervals, J °° dy = dy + dy; in the first interval 
apply the usual approximation 

f dyh(z)f 2 (y) « h(x) f X dyf 2 (y), (C4) 
Jo Jo 

while in the second, since x < y, we can carry out the Taylor expansion as before. Thus, we 
obtain 

I(x) w h(x) + I 2 (x) + h(x) + h(x) , (C5) 

with 

h(x) = A(x) / dyym 2 (y)A(y) , 



o 



j 2 ( x) = / dyi?m*(y)A(y), 



x 

p OG 







/s(x) = -i y ^ 2 A 2 (y)[m 2 (y)]', 
J 4 (x) = -^m 2 (x)x 2 A 2 (x). (C6) 
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Note that, as x — > 0, h(x), hix) and I±(x) vanish, and one recovers from I$(x) the exact 
result for 1(0) given in (IC3I) . 

Finally, note that for a m 2 (x) displaying the asymptotic behavior given in (15.251) the 
following results are useful 



where d is a constant. The first equation is derived using directly the integral of (15.271) . 
while for the other two we have employed the asymptotic property of the incomplete T(a, u) 
function. Specifically, 



(with no restriction on the sign of a), and its asymptotic representation for large values of 
| it | is given by 
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